Persistence in the Zero- Temperature Dynamics of the 
Diluted Ising Ferromagnet in Two Dimensions 
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ABSTRACT 

The non-equilibrium dynamics of the strongly diluted random-bond Ising 
model in two-dimensions (2d) is investigated numerically. The persistence 
probability, P(t), of spins which do not flip by time t is found to decay to 
a non-zero, dilution-dependent, value P(oo). We find that p(t) = P(t) — 
P(oo) decays exponentially to zero at large times. Furthermore, the fraction 
of spins which never flip is a monotonically increasing function over the 
range of bond-dilution considered. Our findings, which are consistent with 
a recent result of Newman and Stein, suggest that persistence in diluted and 
pure systems falls into different classes. Furthermore, its behaviour would 
also appear to depend crucially on the strength of the dilution present. 



In the non- equilibrium dynamics of spin systems at zero-temperature we 
are interested in the fraction of spins, P(t), that persist in the same state 
up to some later time t. For homogeneous ferromagnetic Ising models in 
^-dimensions, P(t), has been found to decay algebraically [1-4] 

P(t)~t~ e{d \ (1) 

for d < 4, where 9(d) is the new non-trivial persistence exponent. 
The presence of a non- vanishing P(t) as t — > oo has been reported in com- 
puter simulations of both the Ising model in higher dimensions (d > 4) [3] 
and the q— state Potts model in 2d for q > 4 [5]; this feature is sometimes 
referred to as 'blocking'. Obviously, if P(oo) > 0, we can reformulate the 
problem by restricting our attention only to those spins that eventually flip. 
Hence, we can consider the behaviour of 

p(t) = P(t) - P(oc). (2) 

Although the numerical simulations of the q— state Potts model mentioned 
above [5] seem to indicate that p(t) also decays algebraically, the evidence 
is by no means conclusive. 

By considering the dynamics of the local order parameter the persistence 
problem can be generalised to non-zero temperatures [6-9]. 

It is only recently [10,11] that attention has turned to the persistence 
problem in systems containing disorder. Numerical simulations of the zero- 
temperature dynamics of the weakly-diluted Ising model in 2d [10] also 



reported that P(oo) > 0. In fact, the study in [10] is consistent with 
the presence of three distinct regimes: an initial short time regime where 
the behaviour is pure-like; an intermediate regime where the persistence 
probability decays logarithmically; and a final long time regime where the 
system 'freezes' and P(t) is effectively constant. 

Very recently, Newman and Stein [11] have argued that the 'blocking' [5] of 
spins in systems with continuous disorder is associated with the fact that 
' every spin flips only finitely many times' 1 . As a consequence, in some simple 
Id models p(t) was found to decay exponentially rather than algebraically 
for large times, namely 

pit) ~ e~ k \ (3) 

where k > 0. In contrast, persistence in the weakly-diluted Ising model 
appears to decay logarithmically in the intermediate regime. (Note that 
[10] examines the behaviour of Pit) and not p{t)). 

Clearly, it is of immense interest to establish whether the presence of 'block- 
ing' in a system necessarily implies exponential decay of the persistence 
probability. Howard [12] has found evidence for exponential decay in cer- 
tain non-disordered models with 'blocking' (2d hexagonal lattices and Bethe 
lattices with z = 3 are discussed in [12]). 

To clarify and further investigate the situation, in this Rapid Commu- 
nication we present the results of computer simulations of an Ising model 
containing strong bond dilution. Here we restrict our attention to zero- 



temperature. 

The Hamiltonian of the model we work with is given by 

H = -J2 JiASj (4) 

where Si = ±1 are Ising spins situated on every site of a square Lx L(= N) 
lattice with periodic boundary conditions and the summation runs over 
all nearest-neighbour pairs only. The quenched ferromagnetic exchange 
interactions are selected from a binary distribution given by 

P(J ij ) = (l-p)8(J ij )+pS(J ij -l) (5) 

where p is the concentration of bonds. 

We obtained data for L = 500 and 750 at zero temperature for a broad 
range of bond-concentrations (0 < p < 0.5) on a suite of Silicon Graphics 
workstations and for L = 1000 on a SGI Origin 2000; as the data for 
the different lattice sizes studied are practically indistinguishable, here we 
simply present the results for the largest lattice simulated. 
We begin each run with a random starting configuration of the spins and 
then update the lattice by first calculating the energy change that would 
result from flipping a spin. The rule we use is: always flip if the energy 
change is negative, never flip if the energy change is positive and flip at 
random if the energy change is zero. 

The number, n(t), of spins which have never flipped until time t is then 
counted. As we are working with strongly diluted lattices, it is necessary 
to monitor the value of n(t) after practically each Monte Carlo step. 



The persistence probability is given by [1] 

P(t) = [< n(t) >}/N (6) 

where < . . . > indicates an average over different initial conditions and 
[. . .] denotes an average over samples i.e. over the bond-dilution. For the 
simulations considered in this work we averaged over at least 100 different 
initial conditions and samples for each run. 

We now discuss our results. To examine the decay of the persistence 
probability, in Fig. 1 we plot lnP(t) versus hit for a wide range of bond 
concentrations, p : 0.1 < p < O.b, for a lattice of size L = 1000. The decay 
of P(t) appears to be non- algebraic before 'freezing' occurs. We see that, 
effectively, P(t) = P(oo) for t > t*(p), where the value t*(p) depends on 
the strength of the dilution. Furthermore, the non-zero value of P(oc) also 
depends on p, with the fraction of non-flipping spins increasing monotoni- 
cally with the bond concentration. The increase in -P(oo) with p can be seen 
more clearly in Fig. 2 where we have plotted some additional data at values 
of the exchange interaction not shown in Fig. 1. The numerical values of 
-P(oo) for the different bond concentrations simulated are also displayed in 
Table 1. 

Obviously, when p = all spins eventually flip as the energy change in 
flipping is always zero. For a value of p ^ 0, there will be regions of the 
lattice containing finite clusters where it will cost energy to flip spins. For 
example, an isolated bond connecting two up spins is just such a stable clus- 
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ter. The occurrence of these clusters increases with the bond concentration 
and hence also does the fraction of spins which never flip. This increases 
smoothly to p = 0.5, the bond percolation threshold, where it appears to 
level off. That is, the maximum value of P(oo) ~ 0.46. Clearly, P(oo) must 
decrease eventually for higher values of p as we know that every spin flips 
infinitely many times for the pure model, p = 1 [1, 11]. 

We now consider the non-algebraic decay of P(t) to P(oo). As discussed 
earlier, it is more convenient to work with p(t) from Eqn. (2). In Fig. 3 we 
replot the data displayed in Fig. 1 as lnp(t) against t. The straight lines 
are linear fits to Eqn. (3) after discarding data for short times. It is evident 
from Fig. 3 that p(t) indeed decays exponentially at large times. Hence, we 
confirm that for the strongly diluted Ising model in 2d persistence decays 
exponentially as predicted by Newman and Stein [11]. This is in marked 
contrast to the behaviour for the pure [1-3] and the weakly diluted models 
[10]. 

To conclude, we have presented new data for the zero-temperature dy- 
namics of the strongly diluted random-bond 2d Ising ferromagnet. This 
system exhibits 'blocking' and we find evidence that p(t) decreases expo- 
nentially for large times. The fraction of spins which never flip increases 
monotonically from zero with increasing bond concentration. Our results 
support the suggestion that the decay of the persistence probability can 
be non-algebraic for certain classes of models. Indeed, for the diluted 2d 
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Ising model the behaviour of p(t) would appear to depend crucially on the 
strength of the dilution. 
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FIGURE CAPTIONS 



Fig. 1 

Log-log plot of P{t) versus t for the bond-diluted 2d Ising model for a range 

of bond concentrations, p; the size of the lattice is 1000 x 1000. 



Fig. 2 

A plot of the fraction of spin which NEVER flip (P(oo)) against the bond 

concentraion p. 



Fig. 3 

Plot of lnp(t) against t for different bond concentrations, p. The straight 
lines are linear fits to the data after discarding the initial short time be- 
haviour. 
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TABLES 



Table 1 



p 


P(oo) 








0.025 


0.0708(1) 


0.050 


0.1340(1) 


0.075 


0.1900(1) 


0.100 


0.2390(1) 


0.125 


0.2813(1) 


0.150 


0.3173(1) 


0.175 


0.3479(2) 


0.200 


0.3732(2) 


0.250 


0.4097(1) 


0.300 


0.4331(2) 


0.350 


0.4453(3) 


0.400 


0.4526(3) 


0.450 


0.4559(2) 


0.500 


0.4576(1) 



The fraction of spins, -P(oo), which never flip at various values of the bond 
concentraion, p. 
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Figure 2 
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